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We consider the ground state of a mixture of two pseudospin- 1 Bose gases with interspecies spin 
exchange in a trapping potential. In the mean field approach, the ground state can be described 
in terms of four wave functions governed by a set of coupled Gross-Pitaevskii-like (GP-like) equa- 
tions, which differ from the usual GP equations in the existence of an interference term due to 
spin-exchange coupling between the two species. Using these GP-like equations, we calculate such 
ground state properties as chemical potentials, density profiles and spin density profiles, which are 
directly observable in experiments. We compare the cases with and without spin exchange. It is 
demonstrated that the spin exchange between the two species lowers the chemical potentials, tends 
to equalize the wave functions of the two pseudospin components of each species, and thus homog- 
enizes the spin density. The novel features of the density and spin density profiles can serve as 
experimental probes of this novel Bose system. 

PACS numbers: 03.75.Mn, 03.75.Gg 



I. INTRODUCTION 



Multicomponent Bose-Einstein condensation (BEC) has been an active subject of research in recent years. People 
have considered BEC of a mixture of two different spinless species [l|-[8| and of spinor gases such as spin-1 [§-[l3[ and 
pseudospin- i gases 

EMI. 

More recently, a mixture of two distinct species of pseudospin- ~ gases with interspecies 
spin exchange was investigated theoretically (l~7l - [23| . In such a mixture, there are N a atoms of species a and Nt atoms 
of species b, while each atom has a pseudospin degree of freedom with basis state a , J,. The particle number of 
each species N a — N a f + N a ^, {a — a,b), is conserved, but the particle number of each pseudospin component of 
each species N a<T , (a =t,|), is not conserved because of the spin-exchange coupling between the two species. Note 
that for a pseudospin-i gas, the total spin of each species is always a constant S a = N a /2, while its z-component is 
S az = (JV - N ai )/2. 

The scattering between any two atoms are now pseudospin-dependent. For the scattering between an atom of 
species a incoming with pseudospin o\ and outgoing with pseudospin 04 and an atom of species /? incoming with 
pseudospin 02 and outgoing with pseudospin 03, the scattering length is denoted ^? C3(T2(Xl - Correspondingly, the 
effective interaction is g°^ a<T1 S(r - r'), where ff"f CT3(T2(Tl = 2 7 r?i 2 ^f (T3(T2(Ti //x a( 3, where ^ Q/3 is the reduced mass of 
the two atoms. For convenience, we define the shorthands = g"" acr for the intraspecies scattering of the same 
pseudospin a, g"g = 2g™~ a for the intraspecies scattering of different pseudospins a ^ a, g^, = g^,^,^ for the 
interspecies scattering without spin exchange, g e = g% b SaS for the interspecies spin-exchange scattering. We shall also 
use the shorthands for ^" acr , for 2£"? crS ., for ^ b a , a , (J . In this paper, all these scattering lengths £'s and 
thus the effective interaction strengths <?'s are considered to be positive. Thus the many-body Hamiltonian density 

is [3 [13 

<T(T 

where U aa (r) is the external trapping potential, i/w = i/' QO -(r) is the Bose field operator for species a with pseudospin 
a (a — a, b, a =f,4,), a and a' may or may not be equal. Expanded in terms of an orthonormal set of single-particle 
orbital basis states {(j> arTii (r)}, 

Vw(r) = E aao-,i0a<T,i(r), (2) 
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where a aa .i is the annihilation operator corresponding to ao -,i(r). In Ji, the first two summations of interaction terms 
are the density-density interactions without spin exchange, as studied in previous models of Bose mixtures, the last 
term is the spin-exchange interaction, which causes spin correlation or entanglement between the two species, and is 
also responsible for the novel features discussed in the present paper. 

Under the usual single orbital-mode approximation, Bose statistics and energetics governs that all atoms of each 
species a and with the pseudospin state a occupy the lowest-energy single-particle orbital mode hereby denoted as 
<^a(r( r )) hence in the expansion ([2]) of the field operator ip a(T (r), one only needs to consider one term a aa -4> aiT (r), where 
a atT is the annihilation operator corresponding to <j> aa (r) . Consequently, in each term of the many-body Hamiltonian 
J d r"H(r), there is an integration of a product of single-particle wave functions, which now becomes an effective 
coefficient. Therefore, under the single orbital-mode approximation, the details of Qcr (r) are not needed in describing 
the many-body ground state in terms of creation and annihilation operators or, equivalently, the collective spin 
operators, although such simplification is lost when one goes beyond single orbital-mode approximation. In a broad 
parameter regime, the two species are quantum entangled in the particle numbers of the two pseudospin states or 
collective spins, that is, the two species do not undergo BEC separately, hence the ground state was dubbed entangled 
BEC. 

However, these four wave functions and the corresponding elementary excitations are important physical properties. 
For a uniform system, cf> aa (r) is simply the constant l/\/F2, where f2 is the volume of the system. For convenience, we 
write the mean field value of the Bose field operator (Vw) as ip a(7 e ila "' , where tp a a > 0. This is the so-called condensate 
wave function. The condensate wave function for pseudospin a component of species a is ip aa — r\ a a \J N® a / where 
N® a is the corresponding particle number in the many-body ground state of the system, determined by the many-body 
Hamiltonian, rj aa = ± is a sign. To minimize the energy, the signs of three components can be chosen to be + while 
that of the other one is chosen to be — . 

In a trapping potential, which is an experimental necessity for BEC of cold atoms, the wave function ip aa (Y) 
is dramatically different from a constant. Moreover, the details of Vw(r) provide experimentally very important 
information, as their modular square is just the particle density, which is directly measurable and is a key observable. 

In this paper, we consider such a pseudospin-i mixture in a trapping potential, and find some interesting properties, 
especially the density and spin density profiles of the four lowest-energy orbital modes {<^w}. There had been many 
calculations on such properties in other types of BEC mixtures [3j, l24| , which demonstrated that a trapping potential 
brings significant features absent in a homogeneous system. 

The GP-like equations can be obtained by using the Euler-Lagrangc equation 

^ 9C =0, (3) 



with C = i X^ctcr ^aa^tipaa — (n) j and then substituting id t as pL aa . One obtains 



r, 2 



2„,. 

aa 



2m a 

+g^\^air\ 2 4'aa- ~ ge^* &s i>aAaS = AWVW ■ (4) 



where g e term is due to the spin-exchange interaction, and is a new feature absent in previous models of Bose 
mixtures. The minus sign comes from the requirement that the phases j a(J of the four components should satisfy 
cos (7^ — 7aj, — 7fc^ + 7/4) = — 1 for the minimization of the energy. For g e > 0, this spin-exchange interaction is like 
an attractive interaction in some way, counteracting the other interaction terms if the latter is repulsive. However, it 
is a new effect, as it depends on the wave functions rather than the densities. 

When g e term is negligible, the system behaves like the usual Bose mixtures with repulsive interactions. When 
g e term is dominant, the system behaves in a way similar to an attractive mixture with intraspecies interac- 
tion negligible. Moreover, to minimize the spin-exchange interaction energy g e $\^] ) ] r '4 } b-x'>Pai + ^i"^^^^^) — 
-Zgei'ati'aiipbiipbt cos(7 at - 7 a ; - 7^ + j H ) = -2g e ^ a ^ a ^ b ^ b -x , i>ati>aiipbiipbf should be maximized, as a kind of 
interference effect, which means that for each species a = a, b, ip a ^-(r) = ip a i(r), hence the density profiles for the two 
pseudospin components of each species tends to be the same. To see this, note that each wave function can be taken 
to be real and positive [25j | . The other terms in the Hamiltonian of course break this equality, as will be studied later 
in this paper. 

Below we shall describe the finding that the larger the interspecies spin exchange g e is, the stronger the overlap 
between the density profiles of the two pseudospin components of each species is. Hence density profiles are very good 
experimental probes of the underlying interspecies correlations. On the other hand, by comparing the experimental 
and theoretical results on the number density and spin density profiles, one may estimate the spin-exchange interaction 
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strength g e . Experimentally, by studying the the effect of g e on the density profiles, one can obtain the information 
on such a mixture. 

Up to now, there is not yet a report on experimental studies of such a spin-exchange mixture between different 
species. However, the interspecies spin-exchange interaction is determined by the difference between the interspecies 
triplet and singlet scattering lengths, which has been found to be quite a few nanometer (nm) [26| . In this paper, the 
theoretical investigation using this parameter value clearly indicates interesting new features. Hence our work also 
provides some motivation and methodology for experimental exploration of such a mixture. 

In Sec. UH the numerical method is described. In Sec. IIII1 we calculate the ground state properties, comparing the 
cases with and without spin exchange, and demonstrating the experimentally observable effects of interspecies spin 
exchange. Then we make a summary in Sec. IIVI 

II. NUMERICAL METHOD 

We assume the trapping potential to be 

U a(T (r) = ^M a cj 2 a (p 2 + X 2 z 2 ), (5) 

where p = yj ' x 1 + y 2 , A represents the trap anisotropy, M a and w Q are the mass of a-atom and trap frequency, 
respectively. U a ^ — U a ±. For a magnetic trap, M a w^ = JoPbBq, with j a being the g factor of a-atom, Bq 
being the central magnetic field multiplied by a normalized factor. In order that the parameter values are close to 
the experimental data, we imagine species a as 87 Rb and species b as 23 Na, then M a oj 2 /Mi,uj 2 = j a /jb = 1, i-e. 
Uaa = Uba- Define k = L)b/ui a y/ M a /Mb = y87/23. In our calculation, we use the parameter values u> a — 2n x 75Hz, 
Ub = K*u> a and A = v8- The values of scattering lengths are set to be ^ = 8 nm, £££ = 6.7 nm, = 3.7 nm, £^ = 4 
nm, £jj = 2 nm, £g = 1.7 nm, ^ = 1.2 nm, = 0.67 nm, £g = 1.9 nm. We vary the value of £ e and study its 
effect on the number densities and the spin densities. 

We expand i/) aa in terms of Nbasis eigenfunctions of the non-interacting Schrodinger equation in an anisotropic 
harmonic potential ([5]), that is, 

N basis 

Vw(r)= Y Ar a eK;{p)®mMZ Wr {z)> (6) 

where R™ r (p), $> mr (tp) and Z Wr (z) correspond to the cylindrical coordinates p, tp and z, respectively, A r aa is the 

expanding coefficient under the condition ^2^ ba " s (A^) 2 — N aa . For the ground state, only eigenfunctions with 
m r = are relevant. 

Therefore, GP-like equation ((4]) is transformed to the following nonlinear equation 
(E l aa - p aa )A l aa + g% A^ aa A k aa I{ia, ]a , ka, la) 

ijk 

Yl ^A^lJiia, ja, ka, la) + gZ Y A L4^L^(^, ja, ka, la) 

ij k ij k 

+9at Y A ^ A iff A aJ(ia, ja, ka, la) - g e Y ^sMta^ff 1 ^ 3", la) = 0, (7) 

ijk ijk 

where E l aa = ((2n ; + m ; + 1) + (wi + l/2)X)huj a , 



1 f°° 

I(ta,ja,k(3,W) = —j JC/(^ Q P 2 )^(^P 2 )<H^P 2 )C'(^P 2 )^P 

/+oo 
Z Wi (j]za z )Z Wj (r] za z)Z Wk (rj z fjz)Z W[ 
-oo 

where m r (r — i, j, k, I) should be 0. There are various algorithms to determine the solutions of the nonlinear equation 
set, such as fixed-point iteration and Newton method. We use Broyden method to obtain the solutions of Eq. ([7]) 
because of its high speed and precision. 
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FIG. 1: Reduced chemical potential p aa /Tiuj a varying with the atom number N of each species, at a generic parameter point 
without interspecies spin exchange, i.e. £ e = 0. p aa as well as paa/^a — Pba/hcdb increase with N. 

ILL CALCULATIONS 

We shall use the GP-like equations (HJ to study the effect of interspecies spin exchange on the number densities 

n a a = |i/w(r)| 2 , 

and the spin densities 

5 QZ (r) = l[|^ T (r)| 2 -|^(r)| 2 ]. 

For simplicity, we assume that the atom numbers of the two species are equal, i.e. N a = N b = N. Although the 
cases of N < 1000 may not be experimentally realistic, they serve as theoretical demonstration of the dilute limit, in 
which the interaction energy is small, especially in comparison with the cases with larger N. 

A. The case without interspecies spin exchange 

First we consider the case of £ e = 0, in which the system reduces to a mixture of four spinless condensates. The 
results for this case are summarized in Figures [H [2J [3J HI El El They are all consequences of minimizing the repulsive 
interactions among the four components, as detailed in the following. 

We have > under the present parameter values, as indicated in Fig. [TJ which shows the chemical potential 
H aa as a function of N. \i aa as well as p a(T /TiLu a — pba/^b increase with N. The difference between and /i a i 
is due to spin dependence of various scattering lengths. This result can be confirmed by a calculation based on the 
Thomas- Fermi approximation, which leads to p aa = U aa dr + g"£N aa + g^N aS + g"^N & a + g^N & <t]- In a 
sense, p aa represents the average energy of one atom of species a with pseudospin a. We can observe that p aa /hui a 
decreases with N, towards the single particle value 1 + A/2 = 2.414. 

Now we take a look at the spatial dependence of the atom densities n a<7 (r) = l^o^r)! 2 . The density profiles 
for several different values of N are shown in Fig. (5J where the distributions along p and z directions are depicted 
respectively. These plots display some interesting features. Obviously n a f and n Q j. are complementary, because of 
the normalization condition f(n a f + n a ^)dr — N . When N is small, n a -\ and n a ^ are close to each other, because 
the interaction energy is small. But with the increase of TV, the difference between n Q -|- and n a x increases in order to 
lower the interaction energy. When N is large enough, two or more peaks may appear in some density profiles, due 
to the inclusion of higher order eigenfunctions in the expansion ((5]) when the interaction energy becomes important. 

The profiles of the total density of each species, n a = n a f + n a ±, are shown in Fig. [3] We can see that not all 
components co-exist in every region, in analogy with the two-component BEC [3]. This is in order to minimize the 
total energy under the given parameter values. 

The density profiles are more extended in p direction than in z direction, as exhibited in Fig. [5] and in Fig. [3] The 
reason is that the trapping in z direction is stronger than that in p direction. 
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FIG. 2: Density n aa = |Vw| 2 for each pseudospin component of each species, at a generic parameter point with £ e = 0, as 
defined in the text, for N = 100, 1000, 5000, 10000. N is atom number of each species. The upper plots are profiles along p 
direction with z = 0. The lower plots are profiles along z direction with p — 0. More than one peak appear in some plots, due 
to the inclusion of higher order eigenfunctions in the expansion © when the interaction energy becomes important. 



We have also plotted the two-dimensional density profiles with p and z as the two coordinates, as shown in Fig. [U 
for A = y/8, and in Fig. [5] for A = 1, both with N = 10000. The complementarity between the two pseudospin 
components of each species is very clear. In the former case, as A ^ 1, the profiles are asymmetric between p and z 
directions. In the latter A = 1, the profiles are symmetric between p and z directions. 

Now we consider the spin density S az (r) = i[|-0 Q ^(r)| 2 — |?/> Q _|_(r)| 2 ]. The spin density profiles are shown in Fig. [5] 
for N — 100, 1000, 5000, 10000. It can be seen that along p direction, the spin density increases from a negative value 
at the center to a positive value at a certain radius, and then gradually decreases to zero. Both the radius with the 
positive maximal spin density and the radius where the spin density becomes zero increase with TV. This feature can 
be understood in terms of the difference between the densities of f and \. atoms shown in Fig. [21 That the spin density 
is negative in inside regimes while positive in outside regime is because we have assumed > consequently 
more f atoms tend to stay in the outside regime where the density is lower because of trapping potential, while more 
J, atoms tend to stay in the inside regime, in order to lower the total energy. Of course the spin density approaches 
zero for large enough value of p, as the densities of f and J. atoms both approach zero. Along z direction with p = 0, 
the spin density is mostly negative, also because more f atoms than \. atoms stay in the larger p regime. This effect 
weakens with the increase of z, because the trapping potential increases, consequently the difference between the 
numbers of I and f atoms decreases. There is only a very small regime where the spin density becomes positive but 
the values are too small to be visible on the plots. This feature is unlike the that of the profiles along p direction, as 
there must be some regime of p with more f atoms. 



B. The case with interspecies spin exchange 

Now we come to the effect of interspecies spin exchange, i.e. the case of £ e ^ 0. We have chosen £ e =0.53 nm, 1.07 
nm, 2.03 nm, 4.27 nm. As stated in the Introduction, the mean-field spin-exchange interaction energy 
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FIG. 3: Profiles of the total density for each species, along p direction with z = 0, and along z direction with 

p — 0, at a generic parameter point with £ e = 0, as defined in the text, for N = 100, 1000, 5000, 10000. Not all components 
co-exist in every region, in order to minimize the total energy under the given parameter values. The density profiles are more 
extended in p direction than in z direction, as the trapping in z direction is stronger than in p direction. 




FIG. 4: Two-dimensional density profiles on a cross section including z-axis, with p and z as the coordinates, in a generic 
parameter point defined in the text, with £ e = 0, A = ^8 and N = 10000. The profiles are asymmetric between p and z 
directions as A / 1. 
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FIG. 5: Two-dimensional density profiles on a cross section including z-axis, with p and z as the coordinates, in a generic 
parameter point defined in the text, with £ e = 0, A = 1 and N = 10000. The profiles are symmetric between p and z directions 
as A = 1. 



is an interference term and acts like an attractive interaction among the four orbital wave functions in some way, but 
it depends on the wave functions rather than the densities. This interference effect is manifested in the number and 
spin density profiles. 

The spin exchange lowers the chemical potential, as evident in Fig. [Jj which is the numerical result of 
function of N for four values of £ e . It can be seen from the plots that the larger £ e is, the lower p aa is. However, in 
all these cases, dfi aa /dN remains positive. 

The interference among the four orbital wave functions is also manifested in the density profiles, as shown in Fig. [8l 
Fig. [9] and Fig. QJJ1 The spin exchange energy (j8]) is minimized when ip a ^ = ip a i for each species a. It is competitive 
with the other interactions. If it dominates the energy, in order to lower the energy, the wave functions of the two 
pseudospin components of each species tend to be close to each other, compared with the case without spin exchange. 
Indeed, this tendency can be observed in Fig. [5] and Fig. [§1 which show the density profiles for each pseudospin 
component of each species, in p direction with z = and in z direction with p — respectively. The larger £ e is, 
the more dominant the spin-exchange interaction is, then the closer tp^ and ip a ± are to each other. When £ e is large 
enough, the overlapping effect becomes very visible. Nevertheless, the overlap is not complete, because of other terms 
in the energy. 

Fig. [10] depicts the profile of the total density of each species, which can be compared with the plot for the same N 
in Fig. [31 it can be seen that the overlap regime of the two species is also enhanced by spin exchange, because of the 
effect of spin exchange term {8]). 

Moreover, we have also calculated the spin density profile, as shown in Fig. [TTJ for N = 10000. Spin density 
Saz = IdV'afl — IVwl 2 ) is proportional to the difference between the densities of the two pseudospin components, so 
it is a quantification of the overlap between the two pseudospin components. Evidently, with the increase of £ e , the 
variation of the spin density with the radius decreases. In other words, both spin density S za of each species and the 
total spin density S z are homogenized by the spin exchange. When £ e is large enough, each spin density tends to 
vanish. The underlying reason is also because ip a ^ and ipai, whose difference gives the spin density of species a, tend 
to be close to each other in order to lower the spin-exchange energy. Consequently the spin density of each species a 
tend to vanish. In both the case without spin exchange and the case with spin exchange, the location on z-axis where 
the spin density becomes zero is much smaller than that along p-direction. Compared with the case without spin 
exchange, another notable feature is that when £ e and N are large enough, in the spin density profile along z direction 
with p = 0, the regime with positive spin density becomes more visible (Fig. Ill[) . This is because the negative sign of 
the spin-exchange interaction counteracts the trapping potential, even though the trapping is stronger in z direction 
than in p direction. 
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FIG. 6: Spin density profiles for N = 100, 1000, 5000, 10000, in a generic parameter point as denned in the main text, with 
£ e = 0. The spin density approaches zero for large enough p or z, as the densities of f and 4- atoms both approach zero when 
the trapping potential is large enough. 



IV. SUMMARY 



This paper concerns the ground state properties of a mixture of two species of pseudospin- ^ Bose gases with inter- 
species spin-exchange interaction in a trapping potential. We have numerically calculated the four orbital condensate 
wave functions, each of which corresponding to a pseudospin component of each species, by using the GP-like equa- 
tions. We set the atom number of each species to be N. Using these wave functions, the number and spin densities 
are obtained. When the spin-exchange scattering length is zero, this mixture reduces to a mixture of the usual type, 
with four components. Various features appear as consequences of minimizing the density-density interaction. For 
example, with the increase of N, the difference between n a ^ and n a ^ for each species a increases. 

If there exists interspecies spin-exchange scattering, novel features absent in the usual mixtures emerge. As the 
spin-exchange interaction is negative as a consequence of minimizing the energy, it acts like an attractive interaction. 
Nevertheless, it depends on the overlap among the four wave functions. It lowers the chemical potentials and make 
the densities of the two pseudospin components of each species tend to be close to each other, and thus the spin 
density tends to be homogenized, and even tends to vanish when the spin-exchange scattering length is so large that 
it dominates over the density-density interaction. 

Therefore as experimentally measurable quantities, the number and spin density profiles of such a mixture with 
interspecies spin exchange are effective probes of the novel many-body ground state of this system. 
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FIG. 7: Reduced chemical potential jj, aa /TiLO a varying with N for four different values of £ e . The chemical potential is lowered 
by the spin-exchange energy, which is attractive. d^aa/dN remains positive. 
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FIG. 9: Density profile n aC r = |Vw| for each pseudospin component a of each species a along z direction on the line p = 
for N = 100, 1000, 5000, 10000 and several values of £ e . Other parameter values are given in the main text. The larger £ e , the 
stronger overlap among the four wave functions. 




FIG. 10: Profiles of the total density n a for each species a along p direction on the plane z = and on the line p — for atom 
number TV = 10000 and several values of £ e . The larger £ e , the closer the profiles of the two species. 
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FIG. 11: Spin density profiles in a generic parameter point as defined in the main text, for N = 10000 and various values of 
£ e - The larger £ e is, the close the spin density is to 0, even at relatively short distances from the origin. 



